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Abstract. 

Nienhuis' truncated 0(n) model gives rise to a model of self-avoiding loops on the 
hexagonal lattice, each loop having a fugacity of n. We study such loops subjected to 
a particular kind of staggered field which for n — > cxo has the geometrical effect of 
breaking the three-phase coexistence, linked to the three-colourability of the lattice faces. 
We show that at T = 0, for w > 1 the model flows to the ferromagnetic Potts model 
with q = states, with an associated fragmentation of the target space of the Coulomb 
gas. For T > 0, there is a competition between T and w which gives rise to multicritical 
versions of the dense and dilute loop universality classes. Via an exact mapping, and 
numerical results, we establish that the latter two critical branches coincide with those 
found earlier in the 0(n) model on the triangular lattice. Using transfer matrix studies, 
we have found the renormalisation group flows in the full phase diagram in the (T, w) 
plane, with fixed n. 

Superposing three copies of such hexagonal-lattice loop models with staggered fields 
produces a variety of one or three-species fully-packed loop models on the triangular 
lattice with certain geometrical constraints, possessing integer central charges < c < 6. 
In particular we show that Benjamini and Schramm's RGB loops have fractal dimension 
Di = 3/2. 



1 Introduction 



The 0(n) model plays an important role in the realm of two-dimensional critical be- 
haviour m 121 El E] • In its lattice version, the partition function at temperature T reads 

Zo(n) = ^ n exp (lf;S, ■ , (1.1) 

where the Si are n-component spins living on the vertices i of some regular two-dimensional 
lattice. The set of nearest-neighbour vertex pairs (lattice edges) is denoted by {ij). 

From the high-temperature expansion, Z is linked to a model of loops, each loop 
having the fugacity n ^J. On the other hand, the 0(n) model is related to standard 0^ 
field theory, where the quartic term encodes the possibility of loop crossings [HE]- In 
two dimensions, this model possesses a second-order phase transition only for \n\ < 2. 
At the critical point, loop crossings are disfavoured, and in particular the limit n — * 
gives access to the scaling properties of self-avoiding loops yj. 

A related model can be obtained by truncating the high-temperature expansion as 
follows [3 Ej 

Zioo,= T.I[(^ + ^S.-Sj), (1.2) 

the model now being defined on a regular hexagonal lattice. The advantage of this 
truncation — and of the particular choice of the lattice — is that to all orders in the high- 
temperature expansion the loops become strictly self-avoiding. More precisely, one has 

^ioop = T-^E^'^^^ (1-3) 
c 

where the summation is over configurations C of self-avoiding, mutually-avoiding loops 
on a hexagonal lattice of size A*" vertices. By construction, the number of loops passing 
through a given vertex is either zero or one. C denotes the number of loops in a given 
configuration, and V the number of vacancies, i.e., lattice vertices not visited by any loop. 

The model defined by Zioop is exactly solvable (in the sense of the Bethe Ansatz) 
along the curves [H E] 

= 2 ± V2^, (1.4) 

for \n\ < 2. Furthermore, it exhibits critical behaviour with algebraically decaying corre- 
lation functions. Its critical exponents along the curve ()1.4j) "^ are known exactly jHEllZj 
and are believed to reproduce those of the generic 0(n) model at its critical point T^. 
We shall refer to this critical phase as that of dilute loops. The curve ()1.4|) ~ is to be 
understood as a line of renormalisation group fixed points that, for each n G (—2,2), 
attract the whole low-temperature region T G (0,Tc). Also along this curve the critical 
exponents are known exactly [HiHl EI- We shall use the epithet dense loops to refer to 
the corresponding critical geometry. 

However, the solution of 2'ioop along ()1.4p ~ does not describe the generic behaviour 
of the 0(n) model in the low-temperature regime 0^21 • The reason is that the model 
defined by Zi^op is unstable towards the inclusion of loop crossings. When included, these 
make the model fiow to the generic symmetry-broken (Goldstone) phase of the model, 
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which can in turn be described by supersymmetric methods GDI • Despite of this 
fact, 2'ioop defines a very interesting model of self-avoiding loops and merits to be studied 
in detail on its own right. 

Note that the sign of T in (jl.4|) is immaterial. This reflects the fact that, with suitable 
periodic boundary conditions, the number of vacancies is necessarily even. Thus, T = 
constitutes another line of fixed points ^21 1^ , henceforth referred to as compact loops. 
This zero-temperature limit is once again exactly solvable [T^ I15j. and for \n\ < 2 it is 
critical with exponents that are known exactly fSlEl- From a geometrical point of view, 
this limit is particularly interesting since only the subset Co C C of loop configurations 
in which the loops are fully packed contributes to Zioop- 

In this paper we shall study the hexagonal-lattice loop model (jl.3j) in the presence 
of a particular kind of staggered field^ w. To define this field, we label the faces of the 
hexagonal lattice by integers k = 1,2,3 in such a way that any two adjacent faces carry 
different labels. Further define Eq as the set of lattice edges whose two adjacent faces 
carry the labels 1 and 3. Given a loop configuration C, let J-' be the number of edges in 
Eq that are covered by a loop segment. The model to be studied is then defined by the 
partition function 

Z = ^n^TV. (1.5) 
c 

Note that taking w ^ 1 does not break those of the lattice symmetries (translations and 
rotations) that respect the above sub-lattice structure. Thus, in the continuum limit we 
can hope to find further conformally invariant critical points for this model. 

Our main objective is to identify the various critical behaviours of the model (jl.Sp . 
and to study its phase diagram and renormalisation group flows with respect to the 
parameters n, T, and w. To this end we employ analytical arguments based on a Coulomb 
gas construction, exact mappings to related lattice models, and numerical transfer matrix 
results. In particular we show that, for \n\ < 2, the model generically exhibits three more 
critical points than the dilute, dense, and compact loops invoked above. The first of these 
is encountered when w; — >■ oo with T = 0, and we shall identify it as a g = state critical 
Potts model. The remaining two critical points are superpositions of a critical Ising model 
(free fermion) and dilute or dense loops, respectively. 

The second part of the paper deals with fully-packed loop (FPL) models on the 
triangular lattice. The exact solutions of the FPL models on the hexagonal US] and 
the square lattices imply that these models have different critical exponents. The 
reason for this non-universality can be traced back to a difference in the dimensionality 
of the target space of the corresponding Coulomb gases. Indeed, while the usual dense 
and dilute loop phases are described by the continuum limit of an interface model with a 
scalar height, the heights of the FPL model on the honeycomb (resp. square) lattice are 
vectors of dimension two (resp. three) ^HIEI- This motivates the study of FPL models 
on other lattices. 

However, it is known numerically [18J that the standard FPL model on the triangular 
lattice^ is in the same universality class as dense polymers [i.e., as the model (jl.3|) along 

^In section 0] we shall give some justification for the nomenclature "staggered field" and discuss its 
relation to a similar construction in the Potts model. 

^This is also true for a class of decorated lattices that interpolate between the square and the triangular 
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the curve ()1.4|) ]. It was pointed out in Ref. [201 that this type of flow from fully-packed 
loops to the low-temperature phase of the model (|1.3p is indeed the generic scenario. 

To avoid this flow, and to get new interesting critical behavior, one may subject the 
triangular-lattice loops to further constraints than just full packing. One such model, 
known as the RGB model, was introduced by Benjamini and Schramm 211 studied 
numerically by Wilson [22] • In this model, the fully-packed loops are simply prohibited 
to take turns through an angle of ±7r/3 at any vertex; each loop carries a trivial fugacity 
of n = 1. However, the critical exponents of this model do not appear to have been 
determined |22i| . 

Here we study several versions of such constrained FPL models on the triangular 
lattice, and we relate their critical exponents to those of the model (jl.5j) . In particular, 
we establish that the fractal dimension of RGB loops is = 3/2. 

The layout of the paper is as follows. In section |21 we introduce the model ()1.5|) and 
relate it to edge (Tate) colourings of the hexagonal lattice. We also show an exact map- 
ping of our three-parameter model to a twelve-parameter loop model on the triangular 
lattice which has been studied earlier 23j. In section |21 we discuss the structure of the 
transfer matrix, which we use to obtain both analytical and numerical results. Section E] 
gives some motivation for the introduction of the staggered field w, by comparing it to 
a similar construction for the Potts model. In section [SI we fix the notation by reviewing 
the Coulomb gas formalism for the FPL model at w = 1 JHl- We then show that for 
w oo the continuum limit of the FPL model coincides with that of a g = state Potts 
model. The same approach also reproduces known results for the finite temperature 0(n) 
model [Ij . In section we discuss the results for the fixed points, critical surfaces, and 
renormalisation group flows in the full phase diagram in the (T, w) plane, at any fixed n. 
In section |7j we derive exact results for RGB loops on the triangular lattice [22] , and for 
some generalisations thereof. We give our conclusions in section [HI 

2 The model and some of its reformulations 

Our loop model is defined on a hexagonal lattice with a set of special lattice edges Eq. 
If the faces of the hexagonal lattice are labeled by the integers k = 1,2,3 as shown in 
Fig. [TJ then Eq is the set of lattice edges between the faces labeled 1 and 3. On such 
a lattice, configurations C of self-avoiding, mutually-avoiding loops are laid down with 
weights given by Eq. (jl.5j) . 

In the present section we review various exact transformations relating this model 
to other discrete lattice models. These transformations become useful in later sections 
when we dress the Coulomb gas of the model ()1.5p and examine its critical properties. 

2.1 Three-colouring model 

When T = 0, the subset Cq of configurations C which carry non-zero weight in Eq. ()1.5|) 
are those in which every vertex is visited by a loop. We shall refer to this as a fully-packed 

lattices |19) . 
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loops (FPL). 

Such loops are related to the three-colouring model introduced by Baxter p!^. In this 
model, the edges of the hexagonal lattice are covered by colors A, B and C, subject to 
the constraint that no two edges of the same colour meet at a vertex. These colouring 
configurations can be brought into contact with FPL configurations as follows ^3 CS] • 

Attribute an auxiliary orientation (clockwise or ant i- clockwise) to each loop, and as- 
sign to each configuration Cq of oriented loops a weight so that the Boltzmann factor of 
the corresponding un-oriented FPL configuration Co is recovered by summing indepen- 
dently over the two possible orientations of each loop. Parametrising the loop weight n 

as^ 

n = e^'^ + e-^'^ = 2 cos(6A), (2.1) 

this is accomplished by assigning a weight e^*^ (resp. e~^*'^) to each clockwise (resp. anti- 
clockwise) oriented loop. To complete the mapping, a three-colouring configuration 
is identified bijectively with a oriented FPL configuration Cg by letting the colour C 
(resp. B) represent a loop segment oriented from a vertex on the even sublattice to a 
vertex on the odd sublattice (resp. from odd to even). [We have here defined the even 
(resp. odd) sublattice as the vertices whose adjacent edges are dual to an up-pointing 
(resp. down-pointing) triangle.] The colour A means that there is no loop segment on 
the corresponding edge. 

In this way, oriented loops become alternating sequences BCBC ■ ■ ■ or CBCB ■ ■ -, 
depending on the orientation. Also note that violations of the FPL constraint (with 
weight T) correspond to defects in which three A colours meet at a vertex. Finally, to 
recover Eq. p.5|) . a weight w is attributed to each Eq edge carrying colour B or C. The 
correspondence between loops and colours is shown in Fig. ^ 

The non-local loop weights can now be turned into local vertex weights, with respect 
to the three-colouring model, by assigning a weight e*'*' (resp. e"*"*") to each vertex where 
an oriented loop turns right (resp. left). The weights associated with the different vertex 
configurations are shown in Fig. |21 Assuming for the moment free boundary conditions, 
every closed loop turns ±6 times more to the right than to the left, whence this re- 
assignment is compatible with Eq. ()2.1|) . 

The three-colouring model (with no defects, i.e., T = 0) can be mapped to a solid-on- 
solid model with a two-component microscopic height z = {zi,Z2) placed at the vertices 
of the triangular lattice dual to the hexagonal lattice. Details of this mapping can be 
found in Ref . jTH] , and are also shown in Fig. |21 When going from one vertex of the dual 
lattice to its neighbour, moving clockwise around an up-pointing (resp. down-pointing) 
triangle, z is incremented (resp. decremented) by the vector A, B, C depending on the 
colour of the direct-lattice edge being crossed. Adopting the normalisation of Ref. |16j . 
one has 

^-(w')' ""-(-^A)- ""A'WsA)' p-^' 

and since A + B + C = the height mapping is well-defined (up to a global shift) for 
each three-colouring configuration. 



•^We are mainly interested in \n\ < 2, in which case A is real For n > 2, A is purely imaginary 
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Figure 1: Loop configuration on the hexagonal lattice. The lattice faces are labeled 
by 1, 2 and 3 as indicated. Among the lattice edges (in grey), those separating 1 and 
3-faces (broken linestyle) define the set Eq (see text). Oriented loops (in black) live on 
lattice edges with colours B or C. Vertices where three A coloured edges meet break the 
fully-packing constraint. 




-iX -iX -iX 
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Figure 2: The microscopic height increments when crossing the coloured edges, here 
shown for an even vertex. We also give the local vertex weights, corresponding to all 
possible left and right turns of directed loop segments (shown as black arrows). 
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(a) 




(b) 
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Figure 3: (a) Allowed vertices (up to rotations and reflections) of the triangular-lattice 
0(n) model. The vertex weights pi (z = 1, 2, . . . , 11) are defined in the figure, (b) As an 
example, it is shown how the summation over the two possible configurations around a 
2-face, given that the six external legs are uncovered, leads to a vertex on the triangular 
lattice with weight pi = + n. 



It is well-known that in the continuum limit this two-component height description 
reduces to a single-component height description when temperature defects are intro- 
duced (T > 0) |3]. Below, we shall show that this is also true — in a slightly more subtle 
way — when one has T = and w > 1. 

2.2 Decimation of 2-faces 

The model ()1.5|1 can be exactly mapped to a more involved loop model on the triangular 
lattice. To see this, consider shrinking the size of all the 2-faces in Fig. ^ to zero (alter- 
natively this means shrinking all edges not in Eq to zero). The result is a loop model on 
the triangular lattice in which loops are still non-crossing, but no longer self-avoiding, 
in the sense that a given vertex can be visited up to three times by the loops. Up to 
rotations and reflections, there are eleven local vertex configurations, as shown in Fig. |21 

The corresponding weights, denoted as pi, • • • , Piij can be related to the three 
parameters T, w and n by summation over the internal structure (the shrunk 2-faces) 
of each triangular-lattice vertex. Note that all edges remaining after this shrinking are 
now of the Eq type. Therefore, the weight w of an occupied edge can be redistributed by 
assigning ^/w to each of the vertices at its extremities. The result is 



where, as usual, each closed loop carries a weight n. 

The model defined by Fig. Q is the triangular-lattice 0{n) model previously studied 
by Knops et al. ^Hj, for general values of the weights pk- The exact parametrisations of 
eight distinct branches of critical behaviour were conjectured as follows. First, both the 



pi = T6 + n p2=w{T* + l) p3 = w{T^ + T) 
P4 = 2wT'^ p^ = p^ = w'^ pq = pr^ = w'^T'^ 
pg = w'^T pio = pii = 



(2.3) 
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0(n) model itself and a certain Potts model on the hexagonal lattice (the latter having 
both edge and face interactions) were reformulated as spin-one 141-vertex models on the 
triangular lattice. The two 141-vertex models were then supposed to intersect exactly at 
the critical branches. The corresponding critical behaviour was then determined, in part 
by numerical simulations on the critical branches, and in part by analogy with an exact 
solution of the square-lattice 0{n) model |24j . 

It is not a priori obvious whether this classification is complete. Actually this is not 
the case: for T = and w = 1, the model ()1.5j) is known to have a compact loop phase 
[TH| ITT)] , which does not form part of the eight branches listed in Ref . [23] ■ 

The mapping (j2.3j) shows that the three-parameter phase space of the model (jl.5j) is 
a subset of the twelve-parameter phase space of the triangular-lattice 0(n) model. In 
particular, in spite of our incomplete knowledge of the latter model, it seems likely that 
some of the critical branches identified in Ref. [23] are good candidates for further critical 
points of the model ()1.5jl . We shall see below that these expectations are indeed fulfilled. 

It is of course unlikely that a putative new critical point of the model (jl.Sp has vertex 
weights (j2.3j) that exactly coincide with one of the critical branches of the triangular- 
lattice 0(n) model. However, it may still map to a point which flows to one of these 
branches under the renormalisation group (RG). Needless to say, it would be a tremen- 
dous task to elucidate the RG flows of the triangular-lattice 0(n) model; for the model 
(II. 5|) this is however feasible (due to the vastly smaller number of parameters) and will 
be done below. 

3 Transfer matrices 

Before moving on to our analytical results for the model ()1.5|) we describe the construction 
of its transfer matrix. Understanding the conservation laws of the transfer matrix will be 
helpful in identifying the operators that are present in the continuum limit. Furthermore, 
diagonalising the transfer matrix numerically will serve as a check — and sometimes as a 
guideline — for the analytical results. We shall also use the numerical determination of 
the central charge and the critical exponents as an aid in dressing the complete phase 
diagram of the model, in section [HI 

3.1 Basis states and algorithmic details 

When dressing the transfer matrix for the model ()1.5j) a natural first question is how 
to deal with the non-local loop weights n. The solution to this problem is to work in a 
space of basis states containing non-local information about how the loop segments in a 
given time-slice were interconnected at former times 03]. The transfer matrix T time 
propagates the system by updating this connectedness information, and at the same time 
builds up the partition function. 

The standard power method provides a convenient way of diagonalising T. Namely, 
when letting T act repeatedly on a suitable reference state, the result will converge to 
the dominant eigenvector, and the rate of growth of the norm of the iterated state will 
yield the largest eigenvalue (which is non-degenerate by the Perron- Frobenius theorem). 
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The reference state is typically taken as one of the basis states. Care must be taken, 
however, since T is often block-diagonal due to the existence of conservation laws (see 
below). In this case, one must ensure that the reference state belongs to the same block 
as the dominant eigenvector. 

The power method also allows for extracting sub-dominant eigenvalues, by iterating 
several states which are kept mutually orthogonal after each iteration. 

One can think of the entries of T as being indexed by a pair of integers (row and 
column number). Relating these integers to the connectedness information which consti- 
tutes the basis states is technically involved [SI], and is once again linked to the issue of 
conservation laws. Fortunately, this complication can be eliminated completely by char- 
acterising each state by a large integer (which is usually easy) and inserting the states 
generated at each iteration into a hash table. 

Finally, the efficiency of the transfer matrix algorithm is optimised by using a sparse 
matrix factorisation scheme, in which T is written as a product of matrices which each 
add a single vertex to the lattice. 

3.2 Conservation laws 

We now examine the conservation laws associated with the transfer matrices for strips 
of finite width L of the lattice model (jl.5|) . 

On the hexagonal lattice there are two natural transfer (or time-like) directions, T\\ 
and 71, which are respectively parallel (7j|) and perpendicular (71) to one third of the 
lattice edges. In Fig. ^ 7[| is vertical and 71 is horizontal. With periodic transverse 
boundary conditions, the minimal number of dangling edges L for these two choices are 
respectively 3fc and 2k (with k being an integer) in the ground state sector. (In other 
words, as a function of L we can expect the various eigenvalues to exhibit "mod 3" 
and "mod 2" fiuctuations respectively. The ground state sector is such that periodic 
transverse boundary conditions respect the sublattice structure of the hexagonal lattice.) 

In the FPL model (with T = 0), for the direction 7[| the number of edges that are 
occupied by loop segments (i.e., the number of B colours plus the number of C colours) 
is conserved ^3]. For entropic reasons, we can expect the largest eigenvalue to reside 
in the sector with the largest number of configurations, which is then the one with loop 
density 2/3. If there are 21 occupied loop segments in a strip of size L = 3k, then the 
number of states in the sector is 



where q = (2/) !/[/!(/ -|- 1)!] are the Catalan numbers. For the direction 71, the conserved 
quantity is Q = [the number of loop segments on even edges] minus [the number of loop 
segments on odd edges]. This can be easily checked by considering the action of the 
transfer matrix on all possible states of an even edge and the adjacent odd edge. For 
Q = 2q, the number of states in the sector labeled by Q is 






(3.2) 
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Namely, out of the available 2k edges, one has to choose / occupied odd edges and / + 2q 
occupied even edges. These can then finally be interconnected in Q+g ways. The largest 
eigenvalue is expected to reside in the Q = sector. 

Taking T > 0, these conservation laws are modified. For instance, in the 7j| direction 
it is now the difference of the up and down oriented vertical loop segments (i.e., the 
number of B colours minus the number of C colours) which is conserved. 



3.3 Series expansions 

We further note that in the large-n limit, for w < 1 and w > 1, the first few terms of the 
series expansion for the free energy per face can be easily written down: 




Comparing these series expansions with our numerically obtained free energies (see below) 
provides a useful check of our algorithms. 



3.4 Eigenvalues and critical exponents 

The central charge c and the scaling dimensions Xi of various correlation functions can 
be extracted from the transfer matrix spectra in a standard way. Let the eigenvalues be 
labeled as Aq > Ai > A2 > . . . (the first few dominant eigenvalues turn out to be real 
and positive). One can obtain c and the Xi from ISU] 

ML) = /(oo) 

and 

ML) - fo{L) - 

Here the free energy per unit area reads /o = —Aln[Ao)/{LM), where M S> 1 is the 
number of rows of the strip; we similarly relate to Aj. The geometrical factor A reads 
A = 2/\/3 for the Ty direction and A = a/3/2 for the T± direction; it ensures that the 
free energy is properly normalised per unit area. 



TTC 



+ o(l- 



27rxj 



o L 



(3.4) 



(3.5) 



3.5 The case T = 0, w ^ 00 

Consider now the transfer direction in the case of T = and w 00. In Fig. |3]we 
show a small portion of the lattice between two time slices; the transfer direction is here 
understood to be vertical. The Eq edges are marked on the figure, and in any allowed 
configuration these must be covered by loop segments. 

We now focus on the two non-ii^o vertical edges in the bottom of Fig. 01 Using the 
fully-packing constraint — and the fact that all Eq edges must be occupied — it is easy to 
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Figure 4: A small portion of the hexagonal lattice between two time slices; the transfer 
direction is vertical. The Eq edges are shown in black. 




Figure 5: The Potts model as a loop model on the square lattice. The circles represent 
the Potts spins, the solid lines are bonds in the clusters, and the broken lines are loop 
segments. 



see that if these two edges are both unoccupied one arrives at a contradiction. The same 
is true if they are both occupied. Therefore, exactly one of them must be occupied. This 
implies that the only allowed sector of the transfer matrix is the one in which exactly 2/3 
of the vertical edges are covered. 

The ground state sector has 3k vertical edges in each time slice. The number of 
covered vertical edges is therefore 2k, i.e., even. 

4 The staggered field 

In the introduction, it was claimed that the parameter w in Eq. (jl.5|) can be considered 
as a staggered field. To make this more clear, we wish to compare it with a similar 
construction in the Potts model. 

One of the nicest reformulations of the square-lattice g-state Potts model IS ClS db 
loop model on the square lattice It is well-known that the Potts model can be 

transformed into a random cluster model, where the summation over Potts spins is turned 
into a summation over bond-percolation clusters [25j; each cluster connects vertices whose 
spins are in the same state. This can be further transformed into a loop model in which 
the loops bounce off the boundaries (both exterior and interior) of the clusters, and of 
the dual clusters. This construction is illustrated in Fig. 
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The partition function is then proportional to 

^Potts = E (4-1) 

loops 

where L is the number of closed loops and B is the number of edges covered by a bond of 
the clusters. The parameter u = [e^ — 1)/ y/Q, where K is the reduced spin-spin coupling 
in the Potts model. According to Fig. the correspondence between the weight u and 
the local behaviour of the loops at the clusters depends on the sublattice. In this sense, 
is a staggered field. Also note that the selfdual point is at m = 1, and only in this case 
is the model exactly solvable |SH]- 

Now consider the limit q —>■ oo, with u = 1. There are just two dominant config- 
urations, each with the same weight, in which loops of length four encircle the even 
(resp. the odd) faces of the lattice. In terms of the spins these are the completely ordered 
(ferromagnetic) and the completely disordered (paramagnetic) states. Moving u away 
from 1 will favour one of these configurations. Thus, the selfdual manifold u = 1 can 
be considered as a two-phase coexistence curve, even for finite q. For g > 4 this phase 
coexistence is first order: taking u through 1 drives the system through a first order 
thermal phase transition with non- vanishing latent heat For g < 4 the coexistence 
becomes second order, and m = 1 is indeed the (ferromagnetic) critical point of the Potts 
model. Even in this latter case, taking u 1 will favour one of the phases and induce 
an RG flow that will take the system to either of the two reference configurations, i.e., 
to non-criticality. 

Our motivation for introducing the model ()1.5|) is taken from the above scenario. The 
parameter w is a staggered field in the above sense. Consider first the case of T = 0. 
When n — i> oo with w = 1, there is a coexistence between three dominant configurations, 
in which loops of length six encircle the faces labeled = 1, 2, 3 respectively. 

In Fig. ini we illustrate a domain wall between faces k = 1 and k = 3, for n ^ 1. The 
states on either side of the domain wall can be thought of as completely ordered phases 
of a three-state Potts model, with spins in the state k being placed at the centres of the 
small loops. Let us write its hamiltonian as follows: 

^3-Potts = -if E - (1 - '^-.,2) . (4.2) 

The Potts spins cxj = 1,2,3 thus live on a triangular lattice, which is however slightly 
shifted across the domain wall. To ensure full packing, a single long loop is necessary 
to separate the phases. The space taken up by the domain wall is dehmited by the two 
black curves. In the ordered phases, per Potts spin there is one loop of weight n and 
three satisfied couplings K. Therefore^ 

Kr^^logn. (4.3) 

Taking w 1 similarly induces a magnetic field H in the three-state Potts model ()4.2|1 . 
with 

if~31ogw. (4.4) 

''Note that when q 00, the loop model based on the g-state Potts model can be similarly related to 
an Ising model. The simple argument above then leads to the relation -ftTising = ^logf/; this was found 
in Ref. ^27j by a more complicated argument, which involved counting the length of the domain wall. 
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Figure 6: The space between the two black curves is filled up by a domain wall (thick 
grey loop segment) that separates two domains of phase 1 (black loops, on the left) and 
phase 3 (grey loops, on the right). Potts spins lie at the centres of the hexagons, i.e., on 
two triangular lattices, which have a slight mismatch. 
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(a) (b) 

Figure 7: (a) The derivative —df/dw at w = 1 is continuous for n = (b) The 

derivative is discontinuous at w = 1 for n = 20. 



For n finite, the curve w = 1 can still be considered as the three-phase coexistence 
curve. Since the correlation length for w = 1 is known to be finite for n > 2 and infinite 
for n < 2 ^1 ^1 Ej, it seems natural to expect that taking w through 1 will induce a 
first-order transition for n > 2 and a second-order transition for n <2. This expectation 
is confirmed in Fig. [7|by numerically evaluating the derivative of the free energy density 
/ with respect to w, using the transfer matrix T\\ for strips of various widths L. 

The staggered field w can be expected to have a richer behaviour than in the Potts 
case. Namely, for w > 1 the three-phase coexistence is broken down to two-phase coex- 
istence (rather than just to a trivial one-phase "coexistence", which happens for w <1). 
This means that it is still possible to find some critical behaviour for i;; > 1; in sectional 
we shall describe exactly how this comes about. 

Consider finally the case of finite temperature T > 0. The parameter T acts so as 
to avoid full packing, while w tries to enforce it. There is thus some hope that the 
competition between w and T may lead to interesting multicritical behaviour. This gives 
rise to a rather rich phase diagram, which will be investigated in section IHl 



5 Coulomb gas formalism 

The critical phases of loop models can be described in terms of effective field theories. 
Exploiting the various symmetries at both the microscopic and coarse-grained levels, one 
can construct such theories, and obtain the central charge and critical exponents exactly 
via the Coulomb gas correspondence [23 UH HI] • 

In section IF!T] we begin by reviewing the case of compact loops (T = and w = 1), 
following Refs. Apart from fixing the notation, this will also serve as a basis for 

section in which we obtain the central charge and the critical exponents for w > 1. 
In section 15.31 we show how the same formalism gives back the known results j3] for the 
case of T > and w = 1. 
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5.1 FPL model at w = 1 



The effective field theory for the long-wavelength behaviour of the model is written in 
terms of a coarse-grained height field, obtained from the microscopic height z already 
defined above. 

A typical configuration of the critical oriented FPL model, considered in terms of the 
colouring model, has domains of ideal states, which are colouring states having the least 
possible variance of z. In other word, ideal states are macroscopically fiat. In an ideal 
state, all loops have length six and the same orientation: as each loop can be oriented in 
two ways, there are six such states. To each of them we assign a coarse-grained height 
h = {hi, /i2) = {z), which is the average microscopic height over a \/3 x -\/3 unit cell of 
the colouring. 

The dominant contributions to the free energy are bounded fluctuations around these 
fiat states, as they maximise the local entropy density. In the continuum limit, the 
height is assumed to be a smooth function of the 2— dimensional coordinates {x^,x'^). 
The part of the partition function accounting for the large scale fluctuations, namely 
Z> = / Phexp(— S'[h]), has an Euclidean action S with three terms: 

S = Se + Sb + Sl. (5.1) 

In Fig. IHl we show the positions of the ideal states in the coarse-grained height space. 
The hexagonal lattice with sides | having the ideal states at the vertices is referred to as 
the ideal state graph I. The nodes of X that correspond to the same ideal state, form a 
triangular lattice with sides 1, and is referred to as the repeat lattice TZ. 

Below we list some important height differences associated with symmetry transfor- 
mations between a pair of ideal states, the order being that of increasing norm. The 
vectors = |(A — C) and m.^ = |(B — A) (the superscript c indicates 'colour') have 
norm |; they are associated with elementary colour exchanges A ^ C and A ^ B, 
respectively; these vectors span a hexagonal lattice. The vectors = C and nig = A 
(the superscript r indicates 'rotation') have norm they are associated with transfor- 
mation of one ideal state to another via 120° spatial rotation; they span a triangular 
lattice. Finally, the vectors mf = (C — B) and = (A — C) (the superscript R 
stands for 'repeat') with norm 1 form the basis of the triangular repeat lattice TZ; they 
are associated with transforming the ideal states back on to themselves. 

Each of the terms in Eq. (jS.lj) can be interpreted from geometrical and symmetry 
considerations. 

The first term is due to the elastic fluctuations of the interface, and has terms 
in gradient-square of the height field. Imposing translational and rotational symmetry 
at the coarse-grained level, the general form of this term must be / d^x {Kii{dhiY + 
-^22(^^2)^ + Ki2{dhi) ■ (9/12)}, where d = {81,82) is the usual gradient. Due to colour 
exchange symmetries at the microscopic level, there are some constraints on these cou- 
pling constants. First, the symmetry B ^ C changing the loop orientation corresponds 
to Z2 ^ —Z2, zi ^ Z\ in Eq. ()2.2j) and keeps the weights in Fig. |21 invariant if we let 
A ^ —A. This implies that Kx2 = 0. Second, the cyclic exchange of colours A —>■ B ^ C 
corresponds to Zi — -^^2, ^2 '^^i ~ |^2 and leaves the weights unchanged. 
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(C,B,A) m' (A,C,B) 




(a) (b) 

Figure 8: (a) The ideal state graph X shows the relative positions of the different 
ideal states in the coarse-grained height space. (6) The labeling of an ideal state — here 
(C, B, A) — is the ordered list of edge colours around an even vertex, with face labels as 
shown (starting from the rightmost edge, and turning clockwise around the vertex). The 
vectors m^, 1112, m^, m^, mf and are defined in the text. 
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This implies Kn = K22', we shall denote this common value by gn henceforth.^ Thus, 

Se = J d^x exp {gn{{dhf + (9/12)'}) • (5.2) 

The mapping of the loop model to an oriented loop model with local complex weights 
as defined in section |21 does not work for loops which experience the boundary. For 
example if the model is defined on a cylinder, an oriented loop winding around the 
cylinder will pick up a total weight of one since it makes an equal number of left and 
right turns. Summing over the two orientations there is a weight of 2 for the loop, 
whereas, by definition of the model, every loop should have an weight n. To correct for 
this, one introduces a background term Sb in the action: 



S'B = -y rf^x(eo•h)7^. (5.3) 

In the above, TZ is the scalar curvature of the space on which the model is defined, and 
eo is the background electric charge, yet to be determined. 

On the cylinder, the scalar curvature is zero, except at the two boundaries at infinity 
along the cylinder axis, so that S'b = exp(27rieo ■ (h(a;"'^, +00) — h(2;\— 00))); we have 
here taken to denote the time-like coordinate. The unique solution that attributes the 
correct phase e^^*'^ to every winding oriented loop reads 

eo = (0, eo) = (^0, = (^0, ^ cos-^ (n/2)) . (5.4) 

The terms in the action described so far constitute a Coulomb gas [SO] with two 
bosonic fields hi and where /12 is coupled to a background charge — 2eo. A crucial 
role is played by the last Liouville term in the action, which reads 

Si^ = j d2xiy[h(x)], (5.5) 

where exp(— [h(x)]) is the scaling limit of the complex weights at any vertex, namely 
exp(±iA) as discussed before. We list below the microscopic weights W with respect to 
colorings around the vertex using the notation of Fig. |H1^6) and Eq. ()5.4|) : 

W{C, A, B) = W{B, C, A) = W{A, B, C) = +i^eo 

W(C, B, A) = W(B, A, C) = W(A, C, B) = -z-cq (5.6) 

6 

The microscopic operator W{x.) is uniform in each ideal state, and is a function of the 
height h G X. Since it is a periodic function of h(x) it can be written as a Fourier series, 

iy[h(x)] = W^eexp(z27re-h(x)). (5.7) 

In the continuum limit, this becomes a sum of vertex operators, and e are the corre- 
sponding electric charges which belong to a lattice 7^|^ C 71*. Here, 71* is the reciprocal 

^We have here sticked to the notation of Ref. ^HI- To compare our formulae with Ref. ^Jji oi^^ needs 
to replace g by K/2'k and e by e/27r. 
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of the repeat lattice 7^, and 7?.^ is the lattice reciprocal to the lattice of periods of Ty(h). 
We now determine which vertex operators in Eq. (j5.7p are the most relevant; these are 
the only ones to be kept in the action. 

In a Coulomb gas, operators are associated with electric and magnetic charges. Elec- 
tric charges e are linked to the periodicity of the height field, and appear in the vertex 
operators (spins waves) exp(i27re-h), as discussed in Eq. 15.71 Magnetic charges m are the 
topological charges of vortex defects in the height. Within the Coulomb gas formalism, 
the scaling dimension of a general operator with electric and magnetic charges e and m 
is computed as [S!] 

^ (e-2eo) 



a; e, m 



9 



29 



H — m ■ m. 

2 



(5. 



The lattice TZ is spanned by the vectors mf = (0,-1) and mf = (-^,|). The 
corresponding reciprocal lattice TZ* is hexagonal, and we denote its six shortest vectors 
by e^. For future reference, we list here these vectors as well as the scaling dimensions 
(found from Eq. (j5.8j) with m = 0) of the corresponding vertex operators: 
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^R 



A. I 

V3 V2' 
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;ef) 
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x[ 


;ef) 


= x( 


;ef 


x[ 


;ef) 


= x{ 





i+2eo 



4 
3 

I' 



29 



(5.9) 



■2eo 



29 



To find the electric vectors appearing in the expansion ()5.7j) . we note that the loop- 
weight operator w{h.) has a higher periodicity than that of m^. Namely, the microscopic 
weights ()5.6|) of the ideal states are invariant with with respect to 120° rotations, which 
are linked to the periodicity of m'". Thus, the lattice TZw determining the periodicity of 
the weights is spanned by the vectors and nig of Fig. |H1 The corresponding reciprocal 
lattice T^l^j/ is hexagonal, and its shortest vectors e'' and the scaling dimensions of the 
corresponding vertex operators are as follows: 





x{e[) 
x(er,) 



4+4eo 
29 



x[e'r^ 
xie": 



xie 



x{el) 

4-4eo 
29 



4+2eo 

23 
4-2eo 

29 



(5.10) 



We see that the most relevant vector is 64, and it is sufficient to keep only the corre- 
sponding vertex operator in Eq. ()5.7|) . 

It has been argued in Refs. ^21 ^| that since the loop weight does not flow under the 
renormalisation group, the corresponding continuum operator iy[h(x)] must be exactly 
marginal. In other words, the screening charge of the Coulomb gas is the electric vector 
corresponding to the most relevant vertex operator in the expansion (j5.7j) . We therefore 
set x{el) = 2, and as a consequence the coupling constant g gets fixed as a function of 
Co (and hence n): 

9 = l-eo. (5.11) 

By virtue of Eq. ()5.1H) . all the critical exponents and the central charge for the system 
now gets determined exactly.^ 

^The compactification radius of the boson /12 is the repeat lattice unit vector length {— 1) and does 
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The central charge of the critical FPL model is obtained by noting that there are two 
bosonic fields hi and /12, where hi is free and /i2 is coupled to the background charge 
— 2eo. This is given as jSHUHl 

6e^ 

c = 2 + 12x(eo,0) = 2 ^. (5.12) 

1 - Co 

The thermal dimension xt, describing the algebraic decay of the energy-energy cor- 
relator, is linked to a magnetic charge. Namely, thermal fiuctuations breaking the FPL 
constraint create pairs of vertices not covered by loops, as shown in Fig. Efa). In the 
height model these are vortices with magnetic charge m = 3A = (v^, 0) and zero electric 
charge. Eqs. ()5.8|) and 1)5.111) immediately give, 

XT = ^il-eo). (5.13) 

Note that the electric charges and e*" that could have been candidates for the thermal 
operator are ruled out: the former because it does not respect the periodicity of the 
weights, and the latter because it is irrelevant (and in particular less relevant than ()5.13|) ). 

Similarly, the one-string and the two-string defects, shown in Fig. M^b) and (c), 
correspond to vortices with magnetic charges mi = 2 A -|- B = (-^,1) = and 
m2 = A -|- 2B = (0, 1) = — respectively.^ At each vortex core there is an addi- 
tional electric charge eg, whose role is to cancel the spurious phase factors arising when 
the strings wind around the core jSHj. The corresponding scaling dimensions xi and X2 
thus read 

2 2(1 -eo)- ^ 

The generalisation to /c-string defects is given in Ref. |16j. Note that the correlation 
function of a two-string defect propagating from to r is proportional to the probability 
of having a loop passing through the points and r, and that the two-string dimension 
is related to the fractal dimension Di of the loop via the relation ES] , 

Di = 2-X2. (5.15) 



5.2 FPL model for w 7^ 1 

Having reviewed the known results for the FPL model with = 1, we now show how 
the results get modified for w ^ 1, primarily because of the change of symmetries. We 
expect that the fixed point aX w = 1 is repulsive, and since the symmetry-breaking 
effect of taking w ^ 1 does not compete with any other parameter (note that we still 
have T = 0) is seems unlikely that the RG fiow will halt at any finite value of w. We 

not appear explicitly in the dimension formula (Eq. H5.8|l ^ due to our choice of the units of the vectors 
A, B and C. But even if such a parameter were retained, the product of it with g would get fixed as in 
Eq. fOT|l . 

^Note that in order for the height mapping to be globally defined, the total magnetic charge of the 
vortices must vanish. In particular, this implies that the vortex and the anti-vortex that participate in a 
two-point correlation function of the thermal or the one-string type must reside at opposite sublattices. 
This is not so for the two-string type correlation function; see Fig.Et^). 
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+3A 




(a) 



-3A 





-(A+2B) 



(b) 





\ 



-(2A+B) 



(c) 



A 

-v^ ■ 

+(2A+B) 



Figure 9: Pairs of thermal defects (a), two-string defects (6), and one-string defects (c) 
at the positions and r. The dashed hnes denote the paths of open BCBC ■ ■ ■ strings. 
In (c) the defect string winds around the vortex core (see text). 
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therefore anticipate a flow towards fixed points at w = and w = oo respectively; these 
expectations are confirmed by our numerical results and the subsequent analysis. 

The fixed point at w = is trivial, as the only allowed FPL configuration is the one in 
which small loops turn around the 1-faces, avoiding all the Eq edges; cf. Fig. Clearly, 
this state cannot be critical. 

The situation for w — cxd is more interesting. Even though the loops have to cover 
all the Eq edges, and stay fully packed, there are many available configurations. We shall 
see below that the ideal state graph X gets fragmented, new electric vectors emerge, and 
that for \n\ < 2 the w = oo fixed point is in the same universality class as the critical 
q = 71^ state Potts model. 

We begin by examining the ideal state graph in Fig. |H1 In the limit w ^ oo, the 
two ideal states (A, 5, C) and (A, C, B) are suppressed, as they correspond to colours 
A lying on top of the Eq edges. We have however checked that for a large lattice with 
suitable periodic boundary conditions there exists a sequence of microscopic updates 
(permutation of two colours around a closed loop) that link the four remaining ideal 
states: 



{C, B, A) ^ (C, A, B) ^ (5, A, C) ^ (5, C, A) ^ {C,B,A) ^ ■ ■ ■ , (5.16) 



where all the intermediate configurations in the sequence have the loops covering all the 
Eq edges. Indeed, transforming (C, A, B) into {B, A, C) requires interchanging colours 
B and C, and this is easily achieved by successively fiipping the direction of all the BC- 
type loops. Transforming (5, A, C) into [B, C, A) is more interesting. Flipping first an 
y4C-type loop of length 6 will create a i?C-type loop of length 18. This latter loop can 
be made bigger by fiipping further AC-type loops that touch its perimeter. Continuing 
in this way, the 5C-type loop will eventually go to the boundary of the system (and thus 
annihilate, due to the periodic boundary conditions), and the region at its interior will 
transform into the ideal state {B, C, A), as required. 

This means on one hand that the two-dimensional ideal state graph in Fig. IHl^a) 
becomes fragmented, and on the other hand that the four remaining possible ideal states 
forming a valid quasi one-dimensional ideal state graph. The repeat lattice is now one- 
dimensional, and is spanned by the vector mf, as shown in Fig. Iiur a). 

Let us consider what this change implies for the three terms in the action (|5.1|) . Of 
the various symmetries used to reduce the number of parameters in the elastic tensor, 
only the cyclic exchange of colours A B ^ C fails. This means that Ku and K22 
can now be different. In fact, the structure of X, as shown in Fig. irUT a). suggests that 
the fluctuations of the flrst height component hi is suppressed, and Ku —>■ 00. Denoting 
K22 = 7Tg the elastic term therefore takes the form 



The argument giving the boundary term Sb is unchanged. So the the bosonic fleld 
h2 is again coupled to the background charge — 2eo given by Eq. ()5.4|) . 

The Liouville term S'l is still given by Eqs. ()5.5|) and ()5.7|) . but the microscopic vertex 
weights ()5.6|) are replaced by: 




(5.17) 



W{C,A,B) 



W{B,C,A) = +z-eo 
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(C,B,A) 




(a) (b) 

Figure 10: (a) The quasi one- dimensional ideal state graph of the T = 0, w ^ oo model. 
(6) The state {A, B,C) is disallowed as no loop passes over the edge shared by the faces 
labeled 1 and 3. 
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W{C,B,A) = W{B,A,C) = -t'^eo 



W{A,B,C) = +z^eo-lnH 
o 

W{A,C,B) = -t'^eo-\n\w] (5.18) 
6 

The crucial thing to note is that the invariance under 120° rotations of the ideal states 
has been lost. 

All of the symmetry-related lattices (the ideal state graph, the repeat lattice, and so 
on) should now be considered as one-dimensional, as should the electric and magnetic 
charges. To see this, consider writing the analogue of Eq. ()5.8p with couplings Kn ^ K22- 
In the limit Kn — >■ cxd, the first component of e = (61,62) becomes immaterial, and any 
non-zero first component in m = (mi,m2) would lead to an infinite scaling dimension, 
whence the corresponding operator is infinitely irrelevant. The change from a two- to 
a one-dimensional lattice will actually modify the computation of electric charges, even 
if in some cases they are still reciprocal to the same height periods as before. To avoid 
confusion, we shall therefore denote the one-dimensional electric charges by the letter G 
in the following. 

The repeat lattice is now spanned by mf = (0, —1). The shortest vectors in TZ* and 
the dimensions of the corresponding vertex operators read, by Eq. (|5.8p : 

Gf = (0,-1) x(Gf) = ifa. 

Gf = (0,l) x(G^) = i^ ^^-^^^ 

The next-shortest vectors in TZ* are G^ = 2Gf = 2(0,-1) and G^ = 2G^ = 2(0,1). 
These are also the shortest vectors in 7?-^. To see this, consider the one-dimensional 
projection of the ideal state graph shown in Fig. ^| Due to the equality of weights in 
the first two hues of Eq. ()5.18p . the weight operator is periodic with half the period of 
the repeat lattice. 

Interestingly, the electric charges G^ and G2 coincide with e^^ and e\ in Eq. ()5.10|) : 
the other e'^'s are now absent. The marginality of the weight operator therefore still leads 
to xiei) = 2, and the same fixation of the coupling, (7 = (1 — 60), as in Eq. (|5.11|) . 

Since there is now only a single bosonic field coupled to the background charge — 26o, 
we get the central charge as 

6=1- (5.20) 
1-60 

We next turn to the thermal exponent. A height defect of magnetic charge 3 A is no 
longer allowed, as this would mean leaving an Eq edge uncovered. This suggests that 
the energy operator is linked to an electric charge instead. We claim further that the 
appropriate vertex operator is expected to have the same periodicity as the vertex weights 
(I5.18|) . since the latter encode the energy of the various microscopic configurations (note 
that these "energies" are purely imaginary in the oriented loop model). The candidate 
electric charges are given by Eq. ()5.19|) . and they also happen to be the shortest vectors 
of TZ* (and hence correspond to the most relevant vertex operators). From this argument 
one should think that xt = x(Gf^); however, our numerical results unambiguously show 
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that the correct identification is 



XT = x(Gf ) = i±^. (5.21) 

It is seen from Eqs. fj5.20|) - ()5.21|) that the critical FPL model at w oo is in the 
same universality class as the q = n"^ Potts model. The Coulomb gas of the Potts 
model is known to exhibit a so-called charge asymmetry fl^ 1^ . which explains the 
seemingly paradoxical identification of xt with x(Gf), rather than with x(G^) whose 
vertex operator is more relevant. 

The two-string defect dimension X2 associated with the magnetic charge mf = (0,-1) 
and the fractal dimension Df are still given by Eqs. ()5.14|) and ()5.15|) . 

String defects with any odd number of strings, however, do not exist at the fixed point 
under consideration. This is not surprising, since the magnetic charges of such defects 
lead to infinite scaling dimensions, as explained above. Note also that in the standard 
formulation of the critical Potts model as a loop model on the square lattice, strings can 
be identified as the boundaries of Fortuin-Kasteleyn clusters jSHl , cf . section HI Thus, in 
terms of the Potts spins an odd number of string is meaningless. 

Yet another argument for the absence of string defects with an odd number of strings 
can be obtained by considering the transfer matrix in the direction 7[|, in the limit 
w — > oo. In section Em we have showed that the conserved number of occupied vertical 
edges in a time-slice must be even. Clearly, this rules out the possibility of having states 
with a single string, and any number of pairwise connected loop segments. 

We conclude this section by presenting some numerical evidence for the correctness 
of the above predictions. In Fig. ^Jwe show the effective central charge Ceg for the T = 
model, as a function of w, for various system sizes L. The data for different values of n 
all agree to show a global maximum at w = 1. By Zamolodchikov's c-theorem |2H1 this 
means a repulsive fixed point, just as expected. The repulsive nature is further confirmed 
by the fact that the peak gets sharper as the system size increases. 

When w < 1, the RG fiow goes all the way down to w = 0, which is the location of 
the trivial fixed point discussed above. Here, c^s = 0, indicating non-critical behaviour. 

When w > 1, c^s rapidly attains a plateau, which extends all the way to w = cxd. 
This indicates an RG fiow towards a fixed point at w = oo. The nature of this fixed 
point is attractive, as witnessed by the fact that the plateau gets broader and fiatter with 
increasing system size. 

Finally note that the values of the central charge at the fixed points w = 1 and w = oo 
are in excellent agreement with the predictions of Eqs. ()5.12|) and ()5.20|) respectively. 

We have also checked the values of the thermal exponent xt at the fixed point w = oo. 
We here coded our transfer matrix so that only configurations in which loops covering 
all the Eq edges were allowed. This led to a reduction in the number of configurations, 
and hence exploration of larger system sizes. We list in Tabled the exponents (obtained 
from Eq. ()3.5|) ) as a function of system size for different n, and compare them with the 
Xt given by Eq. ()5.21|) . The agreement is very good, in particular when allowing for some 
extrapolation of the finite-size data. 
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(c) (d) 

Figure 11: The central charge c as a function of w for different n: (a) n = 1/2, (6) n = 1, 
(c) n = a/3, and (d) n = 2. All the curves are obtained using the two-parameter fits of 
Eq. (jUD- 
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L = 9 


L = 12 
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L = 18 
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0.5592 


0.5482 


0.5439 


0.5417 
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0.5 
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0.7436 


0.7496 


0.8 


V2 


0.9981 


0.9996 


0.9999 


1.0000 


1.0000 


1 


1.0 


1.5286 


1.3716 


1.3141 


1.2898 


1.2772 


1.25 


0.5 


0.9466 


1.3112 


1.3996 


1.4588 


1.4959 


1.5842882 



Table 1: Thermal exponent xt in the limit w ^ oo. 
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5.3 Finite temperature, w = 1 



Let us now show how the same formahsm developed for the FPL models above, can be 
used to rederive the known results for the w = 1 model at finite temperature |3]. The 
crucial difference with the FPL model is that there is now a finite density of the vertex 
configurations shown in Fig. M^a), where three A colours meet. Since these are defects 
of magnetic charge Ah = 3A = (-\/3, 0), this makes the two-dimensional height model 
break down. More precisely, this suggests that hi fluctuations are suppressed, and that 
the model must be described in terms of the /i2 component only. In particular, the ideal 
state graph again becomes quasi one-dimensional. 

Let us first discuss the generic effective field theory, i.e., at an arbitrary but small 
value of T. The elastic term and the background term in the action are given now for the 
bosonic field /i2 by Eqs. ()5.17|) and ()5.3p . The local loop weight operator gets modified 
due to the new weights of uncovered vertices: 

w{C,A,B) = w{B,C,A) = w{A,B,C) = +i^eo 

w{C,B,A) = w{B,A,C) =w{A,C,B) = -i-eo 

6 

w{A,A,A) = -ln[T] (5.22) 

This implies that the coarse-grained loop weight operator (j5.7|) is invariant under all 120° 
rotations. The most relevant electric vector is again G2 = 64 = 2(0, 1). The marginality 
assumption leads to x{el) = 2, implying g = 1 — cq. This leads to the central charge 
being given again by Eq. ()5.20p . 

The thermal exponent is linked to an electric charge in the Coulomb gas, since going 
to a one-dimensional height has projected to zero the magnetic charge formerly linked 
with a temperature-like defect of three A colours. The most relevant vertex operators 
have electric charge Gf and . However, these are ruled out since they do not have the 
full symmetry of the weights ()5.22j) with respect to 120° rotations. The next available 
choices are Gf and Gf , of which the latter is already in use as the screening operator. 
By Eq. ()5.10|) we therefore have 

XT = x(Gf ) = ^ii±^. (5.23) 

1 - Co 

With < Co < 1, corresponding to 2 > n > —2, this is indeed the correct result for 
dense polymers |4j. Note also that xt > 2, so that temperature is indeed irrelevant at 
this fixed point. 

Note that the above argument did not in any way use the exactly known location 
of the critical point (jl.4|) ~. The temperature enters in various non- universal quantities, 
such as the coefficients We in Eq. (j5.7p . It therefore may happen at a particular value 
of T that Wcr = 0. We conjecture that this happens at the critical temperature Tc, 
given by Eq. ()1.4|) "'". The screening charge is then G^, and setting x{Gi) = 2 we obtain 
g = 1 + cq. Hence 

Sen 

c = 1 - — ^, 5.24 
1 + eo 



25 




Figure 12: The schematic phase diagram for the 0{n) model on the hexagonal lattice 
in the {w,T'^) plane at fixed n G (0,2). Grey filled circles indicate critical fixed points; 
for the point at {oo, oo) see the text. The renormalisation group fiows are sketched by 
black arrows. The interpretation of the various critical points is as follows (see text for 
details): C compact loops; D dense loops; d dilute loops; P Potts model; D' tricritical 
dense loops; d' tricritical dilute loops. 



and the thermal operator 

XT = x(G^) = (5.25) 

is now relevant, as it should. Note that Eqs. ()5.24j) - (j5.25|) can also be obtained from 
Eqs. (I5.2()j) and ()5.2Hj) by analytic continuation, i.e., by taking — 1 < eo < 0, which once 
again corresponds to —2 < n < 2 |3]. 



6 Phase diagram 

6.1 Fixed points and renormalisation group flows 

In Fig. El we show the schematic phase diagram for the model ()1.5|1 in the (w,T^) plane 
at a fixed value of n. This diagram is valid for < n < 2. The case n = 2 is slightly 
special and will be discussed separately in section 16.41 

The critical fixed points C (compact loops), D (dense loops) and d (dilute loops), were 
previously known and have been discussed in the introduction. The direction of the RG 
fiows for w = 1 follow from the relevance (resp. irrelevance) of xt at C and d (resp. at D). 
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Note also the agreement with the c-theorem j2Hl; according to which the effective central 
charge decreases along RG ffows. 

In section we have demonstrated that for T = there is a flow from C to another 
fixed point P (Potts model) situated at w = oo. 

It should be noted that the points D and P are both generic in the sense that their 
Coulomb gas is one-dimensional with no tunable parameters. (In particular, their central 
charges coincide.) This would seem to indicate that these points should be attractive in 
all directions. In the case of D this is supported by the irrelevance of the temperature 
perturbation, i.e., > 2 in Eq. (j5.23|) . In the case of P, we have seen in section 
that 1/w is irrelevant. It is also true that T is irrelevant. This can be seen, on one hand, 
microscopically by noticing that the w —>■ oo limit enforces full packing for any finite T. 
On the other hand, one has formally x(0,3A) = oo in Eq. ()5.8|) . cf. the remark after 
Eq. ()5.18|) . Thus finite T is indeed infinitely irrelevant at P.^ 

Accepting the attractive nature of D and P in all directions in the [w, T^) plane 
as a working hypothesis, there must be a curve emanating from C that separates their 
respective basins of attraction. For the qualitative reasons given towards the end of 
sectional it seems reasonable to suppose there will be an RG flow along this curve, going 
from C towards a new multicritical point D' situated at finite T and w. 

The simplest way of obtaining a phase diagram that is consistent with the properties 
listed above is to suppose the existence of yet another fixed point d' at infinite T and w, 
which is repulsive in all directions in the {w, T^) plane. This yields then the conjectured 
phase diagram shown in Fig. ^1 The limit w,T oo should be taken with A = T^/w 
fixed, since uncovering an Eq edge will necessarily lead to the formation of two empty 
vertices. Taking A — ^ we recover an FPL model, i.e., the point P. On the other hand, 
A — i> cx) yields a completely empty configuration. The hypothesis is thus that there exists 
a critical fixed point d' situated at some finite value A = Ac. 

6.2 Universality classes and numerical results 

To verify the hypotheses made above — and to establish the universality classes of the 
proposed fixed points D' and d' — we proceed to make some numerical checks. 

To get some guidance, we show in Fig.[T31a plot of the effective central charge (i.e., 
the value of c obtained by fitting the numerically obtained free energies to Eq. (j3.4j) for 
two different system sizes L) as a function of T, at fixed n = y/S and w = 1.60. At T = 0, 
with increasing system size L we have Ccs 4/5, in agreement with the prediction for 
the point P, cf. Eq. (I5.20j) . At T = T^{w) ~ 0.39, one observes a maximum in Ces- This 
is the signature of the line of RG flows that separates the basins of attraction of D and P 
in Fig. El This is followed by a plateau in Ces signaling the basin of attraction of D. The 
broadening of this plateau with increasing L is consistent with the attractive nature of 

^This statement is seemingly contradicted by the fact that < 2 at P [see Eq. H5.21|l ]. However, we 
remind that the notation xt refers to the scahng dimension of the energy operator, which is related to the 
next-largest translationally invariant eigenvalue of the transfer matrix. This does not necessarily coincide 
with the operator that couples to the temperature variable T. Thus, the exponent xt of Eq. (|5.21|) must 
correspond to one of the nine independent way of leaving the manifold H2.3f) which embeds the model 
(|1.5|l within the more general triangular-lattice 0(7i) model. 
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Figure 13: Plot of Ces versus T at n = ^/3 and w = 1.60, for various system sizes L. The 
two maxima are at T~{w) ^ 0.39 and T^{w) ^ 1.84. 



D. Moreover, the plateau value Ces — 4/5 is in agreement with the prediction for D. At 
T = T~^{w) ~ 1.84, the plateau terminates in a second maximum. This is the signature of 
the line of RG flows that separates the basins of attraction of D and the trivial fixed point 
at (w,T^) = (0, oo). For T > T+, c^s levels off to zero, signaling non-critical behaviour. 
Note that the sharpening of the two maxima at T^{w) with increasing L is consistent 
with the corresponding RG flows' being repulsive in the temperature direction. 

The curves Ces (T) have a similar functional form (two maxima separated by a plateau) 
for other values of w. This enables us to determine the curves T^{w) for any fixed value 
of n. In Table El we show the results for T~(w) at n = ^/S. Extrapolations to the limits 
1(7 —> 1 and w oo confirms that this RG separatrix connects the points C and d' in 
Fig. 1121 Actually, we expect T~(w) to diverge as w^^'^ when w ^ oo, since d' is located 



w 


c(3,6) 


c(6,9) 


c(9, 12) 


T-iw) 


1.60 


1.1541 


1.2660 


1.2895 


0.39 ±0.01 


2.00 


1.1847 


1.2762 


1.2921 


0.56 ±0.01 


3.30 


1.2206 


1.2838 


1.2943 


0.99 ±0.01 


3.42 


1.2204 


1.2839 


1.2943 


1.028 ±0.001 


3.50 


1.2200 


1.2838 


1.2942 


1.05 ±0.01 


4.00 


1.2158 


1.2832 


1.2941 


1.19 ±0.01 


6.00 


1.1919 


1.2786 


1.2930 


1.66 ±0.01 



Table 2: Low-temperature maxima of Ceff(T), as a function of system size L, for various 
values of w and fixed n = We also show the estimates of T~(w), extrapolated to the 
L — s> oo limit. 
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at finite A = T"^ /w. Examining now the variation of Cefr (T (w)), we observe a maximum 
at 

Wc = 3.42 ±0.01 and T"(«;c) = 1.028 ± 0.001 (n = Vs). (6.1) 

Tliis is our estimate of tlie location of the point D', here for n = a/3. Note that the 
maximum gets broader with increasing L, in agreement with the attractive nature of 
D' along the separatrix; see Fig. ^| Extrapolating the corresponding central charge to 
L — s> oo we obtain 

CD, (n = ^3) = 1.304 ±0.005. (6.2) 

We have repeated this analysis for n = ^/2, finding similar qualitative results. The 
point D' is here located at 

= 3.60 ± 0.01 and r-(«;c) = 1.030 ± 0.001 (n = ^2), (6.3) 

and the central charge reads 

CD, (n = v/2) = 1.004 ±0.005. (6.4) 

The mapping from the model (jl.5|) to the triangular-lattice 0(n) model j2Sl exhibited 
in section 17!^ makes it plausible that the points D' and d' should coincide with critical 
behaviour found in the latter model. We therefore conjecture that the universality class 
of D' is that of a superposition of dense loops and a critical Ising model, referred to as 
branch 5 in Ref. [23 ■ This means 

3 6el 
2 ~ 1 - eo 

XT = 1, 

since the thermal exponent of the Ising model is always more relevant than that of dense 
loops. 

This conjecture for c is in excellent agreement with Eqs. (|6.2p and (j6.4|) . The numerical 
estimates for xt{L) read, for n = ^/S, 

xt{6) = 0.8622, a;r(9) = 0.8782, xt(12) = 0.8880. (6.7) 

For n = \/2 we find 

0:^(6) = 0.9920, xt(9) = 0.9973, xt(12) = 0.9986. (6.8) 

The agreement with Eq. ()6.6|1 is reasonable for n = ^/S (the finite-size corrections are 
strong), and excellent for n = \/2. 

Following the second maximum in Fig. El allows us similarly to trace out the curve 
T^{w). The result is shown in Table El This separatrix is seen to connect the points d 
and, quite plausibly, d'. Once again we expect T~^{w) to behave as w^^'^ when w ^ oo, 
and this is confirmed by the data in Table El Furthermore, the variation of Ces with L is 
now consistent with d being attractive along the separatrix, exactly as shown in Fig. El 

We have also followed the maxima in Fig. El into the region < w < 1. We find 
that for a certain value w = wq the two maxima coalesce, and for w < wq the model is 
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(6.5) 
(6.6) 



w 


c(3,6) 


c(6,9) 


c(9, 12) 


T+{w) 


1.04 


0.859630 


0.856344 


0.856625 


1.6079 ±0.0001 


1.08 


0.859414 


0.856323 


0.856622 


1.6277 ±0.0001 


1.12 


0.859065 


0.856291 


0.856616 


1.6471 ±0.0001 


1.20 


0.858011 


0.856190 


0.856599 


1.6847 ±0.0001 


1.60 


0.847764 


0.855149 


0.856392 


1.845 ±0.001 


2.00 


0.833671 


0.853567 


0.856042 


1.967 ±0.001 


4.00 




0.843420 


0.853420 


2.456 ±0.005 


6.00 




0.836064 


0.850333 


2.80 ±0.02 



Table 3: High-temperature maxima of Ccfr(T), as a function of system size L, for various 
values of w and fixed n = a/3. We also show the estimates of T~^{w), extrapolated to the 
L — > oo limit. 



n 


c(3,6) 


c(6,9) 


c(9, 12) 


c(12,15) 


Ac 


2 


1.304 


1.409 


1.411 




0.79 ±0.01 




1.222 


1.344 


1.367 


1.369 


0.83 ±0.01 


V2 


1.084 


1.207 


1.241 


1.253 


0.88 ±0.01 


1 


0.825 


0.928 


0.961 


0.976 


1.00 ±0.01 



Table 4: Maxima of Ceff(A) in the limit w,T ^ oo with fixed A = T'^/w, as a function of 
system size L, for various values of n. We also show the estimates of the critical point 
Ac, extrapolated to the L — oo limit. 

non-critical for any T. The interpretation in terms of RG flows is shown in Fig. El Note 
in particular that it is possible to flow directly from C to d, by perturbing C by a flne- 
tuned linear combination of the operators w and T^, and that the dense-loop phase (i.e., 
the basin of attraction of D) only exists for w > Wq. For n = v^, we have determined 
wq = 0.32 ± 0.01, corresponding to a temperature Tq = 0.92 ± 0.01. 

Let us flnally examine the properties of the point d'. As discussed above, we must 
take the limit w,T ^ oo with A = T'^/w flxed in order to see competition between FPL 
conflgurations (with weight ~ w'^^'^, where N is the number of vertices of the lattice) 
and empty space (with weight T^). In this limit, the participating conflgurations consist 
of loops along which Eq and non-i?o edges alternate. An EQ-edge not covered by a loop 
carries a weight A. This reformulation allowed us to write a transfer matrix for this limit, 
with parameters A and n (the loop weight). 

In Fig. El we display a plot of Ceff(A) at n = \/2. As A ^ we recover the point P 
with c = 1/2, and in the large-A limit non-critical behaviour is signaled by Ces = 0. A 
maximum at A = Ac — 0.88 corresponds to the critical behaviour of the point d'. The 
sharpening of the maximum with increasing L shows its repulsiveness in the variable 
A -Ac. 

In Table 01 we locate the point d' for several values of n and we give the corresponding 
values of the central charge. The latter prompt us to conjecture that the universality 
class of d' is that of a superposition of dilute loops and a critical Ising model, referred to 
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Figure 14: Ces as a function of A = T'^/w with T, w — >• oo and fixed n = \/2. 



n 




xt(6) 


xt(9) 


xt{12) 


2 


0.7106 


0.4975 


0.4200 


0.3685 


Vs 


0.8509 


0.6094 


0.5411 


0.4998 




1.1245 


0.7805 


0.7203 


0.6886 


1 




1.1082 


1.0435 


1.0237 



Table 5: Thermal exponent Xt at the point d', as a function of system size L, for various 
values of n. 

as branch 4 in Ref. [22] • In other words, 

_ 3 _ 6el 
^~ 2 " 1 + eo 

We have verified that the point d' can be detected as in Fig. ^^for < n < 2. 

As to the thermal exponent, one would expect xt to be the least of the values given 
by Eqs. ()5.25|) and ()6.6p . Curiously, this is not brought out by the data of Table |31 And 
the numerical data reported in Table 3 of Ref. i23ij for branch 4 of the triangular-lattice 
0(n) model agree neither with the theoretical expectation, nor with the data of Table El 

6.3 Coulomb gas arguments for multicritical behaviour 

Among the various critical branches identified for the 0(n) models on the square (21] 
and the triangular [22] lattices, and for the hexagonal-lattice 0(n) model in a staggered 
field [here], four appear to be generic. Those are dense (D) and dilute (d) loops, and 
their respective superpositions with a critical Ising model (D' resp. d'). 
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(6.9) 



In Ref. inn]; it was explained how these four branches may be recovered in a more 
general setting from Coulomb gas arguments. The reasoning given in Ref. [321 applies 
to the square lattice, but can be readily adapted to the triangular-lattice 0(n) model 
defined in Fig. El 

Indeed, consider a generalisation of the triangular-lattice 0(n) model in which loops 
can be either black or white, with respective weights rzb and n„. Each edge is covered 
by either a black or a white loop. At a vertex, two loops can cross only if they have 
different colours. This can again be illustrated by Fig. IHl but the vertices marked p2 
and p4 now each correspond to two distinct vertices depending on how the white loops 
are interconnected. Clearly, the original model is recovered by setting = 1 and 
renormalising the weights of the "doubled" vertices. 

The Coulomb gas is once again based on an interface model on the dual lattice. After 
orienting each loop, two-dimensional microscopic heights are defined, depending on the 
orientation and colour of each edge. On the square lattice, this construction was based 
on an edge-colouring model [HS] (like in section ^ which defined loops as alternating 
sequences of two colours. This is not possible on the triangular lattice, as the lattice 
is not bipartite and loop lengths are not necessarily even (2D|. The resulting height 
model, however, still exists and has the crucial property that the two loop colours define 
independent directions in the height space^ (this construction — as well as the detailed 
arguments — is exactly the same as in Ref. [SH])- 

The arguments given in Ref. jHS] thus go through, and show that the two loop colours 
lead to two independent one-component Coulomb gases in the continuum limit. There- 
fore, four critical branches arise depending on whether each loop colour is in the dense 
or in the dilute phase. In particular, setting n„ = 1, the white loops add either or 1/2 
to the central charge, depending on whether they are in the dense or the dilute phase. 
Indeed, the dilute 0(n„ = 1) model is nothing but a critical Ising model. 

6.4 The case of n = 2 

The loop weight = 2 is somewhat special. At = 1, the critical points D and d become 
identical and describe the Kosterlitz-Thouless point of the XY model. In the usual 0(2) 
model, this point is the critical end-point of a line of low-temperature fixed points along 
which the exponents vary continuously. 

Interestingly, this feature is not present in the loop model we recall that the 
latter is defined by a truncation of the full 0(n) model, which is not innocuous at low 
temperature. Thus, in the whole region < T < T^{w = 1), the critical exponents 
are constant and equal to their values at T^{w = 1).^° In particular, the thermal 
exponent is marginally irrelevant (resp. marginally relevant) for temperatures slightly 
below (resp. above) T'^^w = 1). 

For w > 1, the lines T~{w) and T~^{w) are present, just as for other n < 2. In 
particular we find the D' critical point at Wc = 2.9 ±0.1 and T~{wc) = 0.91 ± 0.04. Its 
central charge agrees well with the prediction of Eq. (j6.5j) . which reads c = 3/2. 

^This is in contrast with the appearance of an "unexpected" extra height component in the FPL 
models on the hexagonal ^Hj and square ^] lattices, which are both linked to edge-colouring models. 
^"Notice that in this respect the phase diagram reproduced in Ref. 0J is slightly misleading. 
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w 


c(6,9) 


c(9, 12) 


a;i(9) 


xi(12) 


xr(9) 


xt(12) 


T-{w) 


1.6 


1.46271 


1.48616 


0.2277 


0.2249 


0.7151 


0.7231 


0.43 ±0.01 


2.0 


1.47149 


1.48938 


0.2303 


0.2282 


0.7226 


0.7283 


0.60 ±0.01 


2.9 


1.47670 


1.49182 


0.2401 


0.2411 


0.7136 


0.7150 


0.91 ±0.01 


4.0 


1.47452 


1.49078 


0.2401 


0.2407 


0.7051 


0.7067 


1.23 ±0.01 


6.0 


1.46613 


1.48739 


0.2434 


0.2441 


0.6782 


0.6786 


1.70 ±0.01 


10.0 


1.45331 


1.48173 


0.2645 


0.2703 


0.6186 


0.6107 


2.40 ±0.01 


40.0 


1.43237 


1.46268 


0.3106 


0.3169 


0.4967 


0.4655 


5.46 ±0.01 



Table 6: Estimates of Cefr, xi and xt along the T (w) curve at ri = 2, for various system 
sizes. 

Moreover, from Eq. ()6.9|) we also expect a value of c = 3/2 at the point d'; the 
numerical data in Table |3] are consistent with this expectation. This opens for the 
interesting possibility that the part of the curve T~{w) that connects D' and d' may be 
a line of fixed points with c = 3/2. 

The numerical data for c, xt and xi shown in Table IHl are not inconsistent with this 
possibility. While c and xi are fairly constant along the proposed line of fixed points, xt 
steadily decreases beyond w > Wc (also note that the w ^ oo limit of TablelHlis consistent 
with the value of xt given in Table Ej). However, if this scenario were indeed true, one 
should be able to identify an exactly marginal operator at D'. We plan to investigate this 
issue further in the future. 

6.5 Hard hexagons and hard triangles: n > 2 

When w = 1, taking the limit T, n — >■ cxd with fixed z = n/T^ transforms the model (jl.5|) 
into Baxter's model of hard hexagons jlOl EH] ■ The role of the hexagons is played by the 
loops of length six, which may surround any of the hexagonal-lattice faces. As a function 
of their fugacity z there is a solid-to-liquid transition dX z = z^ = (11 ±5^5) /2 at which 
the model is exactly solvable. For ^ ^ oo there is a coexistence of three completely 
ordered phases, in which the loops are fully packed and all surround faces with the same 
label. This may lead one to expect some connection with a three-state Potts model, 
along the lines of section HJ And actually the exact solution dX z = z^. shows |40J that 
the critical exponents of the hard hexagon model coincide with those of the three-state 
Potts model. 

More recently, it was shown numerically jH] that this critical behaviour extends to a 
critical curve Zc{n) existing for any n > 2. In particular, as n — >■ 2_|_ one has Zc{n) — >■ oo 
(or T ^ 0). This remarkable finding means that the symmetry is felt even by the 
"soft" particles (loops of length longer than six) at finite n. The results of Ref. [JIl can 
be interpreted in RG terms by assuming that there is a curve of RG flows along the curve 
Zc{n) directed towards the critical fixed point aX n = oo. Note that an RG flow changing 
the symmetry related parameter n is possible exactly because at n > 2 the loop weight 
operator is no longer marginal. 

In section m we have discussed how taking w > 1 breaks three-phase coexistence down 
to two- phase coexistence. One is thus led to conjecture that for n > 2 there will exist a 
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y\3/\2/\ 

Figure 15: Triangular lattice obtained by superposing three shifted hexagonal lattices. 
The vertices (here surrounded by circles for clarity) and edges are labeled by integers 
k = 1, 2, 3 as shown. 



critical manifold in the w > 1 half space that belongs to the Ising universality class {Z2 
type symmetry). 

Let us consider more closely taking first w ^ 00, and then T,n ^ 00 with fixed 
z = n/T^. Decimating the 2-faces as in section IT^ we recover a model of hard triangles. 
These triangles are loops of length three that surround the faces of the triangular lattice; 
they may touch at the vertices, but not share an edge. The fugacity of a triangle is z. 
The two completely ordered states correspond to having all the triangles surround the 
faces labeled 1 or 3, respectively. It seems natural to conjecture that coexistence of these 
two states may lead to critical behaviour of the Ising type.^^ 

We have checked these conjectures numerically. For any n > 2 we find that at w = 1 
there exists a critical temperature T(.{w = 1) with a c = 4/5 behaviour, consistent with 
the findings of Ref. |^. For w > 1, this extends into a critical curve Tc{w), with an 
RG flow directed towards larger w that eventually ends up in a c = 1/2 critical fixed 
point in the w ^ 00 limit. The location of this fixed point is such that A = T'^/w tends 
to a constant. This flow from c = 4/5toc=l/2 divides the phase diagram into a 
low-temperature ordered (solid) state and a high-temperature disordered (liquid) state; 
there appears to be no further critical points. 



7 Loops on the triangular lattice 

We finally consider several non-standard FPL models on the triangular lattice. These 
loop models exhibit different critical behaviour than the standard FPL model on the 
triangular lattice, which is in the universality class of dense loops [TH} l2Uj. 

Consider constructing a triangular lattice by superposing three shifted hexagonal 

^^In order to compare with the conventional formulation of models of hard particles, it is useful to 
consider the dual lattice. The particles of fugacity z now occupy the vertices of the hexagonal lattice, 
and the hard-core constraint is that no two neighbouring vertices can be occupied. In this sense, the 
model is one of hard squares whose centres live on the hexagonal lattice. Note that this is different from 
the hard squares model where the particles live on the square lattice. Neither of these models appear to 
have been exactly solved. 
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(a) (b) 

Figure 16: (a) The thick grey loops are the RGB loops. (6) The three hexagonal lattices 
are shown by the continuous, broken, and dotted lines respectively. Following successive 
dimers along edges on these three lattices, as shown by the arrows, one gets an RGB 
loop. 



lattices, as shown in Fig. ^1 The vertices of the triangular lattice are labeled k = 1,2,3 
so that no two adjacent vertices have the same label. Let (/ci, ^2, k^) be any permutation 
of (1,2,3). Then the /ci'th hexagonal lattice consist of those edges of the triangular 
lattice that join vertices labeled k2 and ^3. Conversely, at a triangular-lattice vertex 
labeled fci, the hexagonal lattices labeled ^2 and k^ share a common vertex. Each edge 
of the triangular lattice belongs to exactly one of the three hexagonal lattices. 



7.1 The RGB model 

The first model we consider is known as the red-green-blue (RGB) model, and was in- 
troduced and studied in Refs. j2Il 122] • It is obtained by superposing independent dimer 
coverings of each of the three hexagonal lattices. We shall here study the more gen- 
eral case where the hexagonal lattices are endowed with three-colouring configurations, 
cf. section 12. H each BC-type loop being weighted by n. Dimers correspond to colours 
A, and to weigh each dimer configuration equally we simply set n = 1. The dimers now 
form closed loops — the RGB loops — which are each weighted by unity. Apart from being 
fully packed, the RGB loops have the interesting properties that their length is always a 
multiple of three, and that they never make n/S turns. An example is shown by the thick 
grey loops in Fig. ITHT a) (disregard for the moment the other two kinds of loops shown). 

In Ref. 1221 ^ one- dimensional height model was associated with the RGB loops as 
follows. Imagine starting at a definite vertex of the triangular lattice, and traversing the 
RGB loop that goes through this vertex by moving in the clockwise direction. Each edge 
on the trajectory is a dimer on one of the three hexagonal lattices, and when tracing out 
the loop we note sequentially the labels of the hexagonal lattices to which each of these 
dimers belongs. If this sequence of labels is 123123 ■ ■ ■, or a cyclic permutation thereof, we 
shall call the RGB loop clockwise oriented; otherwise it is anticlockwise oriented. These 
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orientations can be used to define a one- dimensional height on the lattice faces, which 
we shall refer to as the RGB height. Properties of this height were studied in Ref. [221 • 

Note however that the orientation constructed above is not an intrinsic property of 
the RGB loops: unlike the orientation of BC-tjpe loops defined in section ITTl the RGB 
orientation is fixed once the trajectory of the corresponding loop has been fixed. For this 
reason we shall in the following rather focus on the orientations of the BC-type loops, 
and on the associated height mapping. 

Since the three copies of the fully-packed hexagonal-lattice 0{n = 1) model do not 
interact, the critical properties of RGB loops are easily found. The free energies and the 
central charges of the three copies simply add up, and in particular one has 

c = 3. (7.1) 



The same is true for critical exponents linked to vortex-type defects. Correlation func- 
tions in the RGB model are products of those in the individual hexagonal-lattice models. 
So once a given defect configuration on the triangular lattice has been decomposed into 
defects in the three hexagonal-lattice models, the corresponding critical exponents just 
need to be summed up. 

In Fig. El we show a few examples. We first consider the thermal exponent on the 
triangular lattice, which is associated with an empty vertex. The corresponding vortex 
configuration, shown in Fig. ITTT a). is equivalent to two one-string defects on the two 
hexagonal lattices which meet at the centre of the vortex. Thus, from Eq. (|5.14j) . 

x^" = 4- + 4'=^ = 1. (7.2) 

Note that this is different from xt of the n = 1 FPL models on the hexagonal and the 
square lattices (which have xt = 1 and xt = 2/3 respectively). Also, in order for the 
energy-energy correlation function to have zero total magnetic charge on each hexagonal 
lattice, the two thermal defects have to placed on triangular-lattice vertices carrying the 
same label. 

Next, consider the bend defect associated with a loop turning by an otherwise forbid- 
den angle of ±7r/3. As all non-defect turns go through angles of ±27r/3, it is easy to 
show that a loop cannot close if there is only one ±7r/3 defect. We henceforth consider 
an RGB loop with two such defects, situated at distant points, and corresponding to 
angles +n/3 and — 7r/3 of opposite sign. The vortex due to one of these defects is shown 
in Fig. ITTT fe): it amounts to a one-string defect on either of two hexagonal lattices, just 
as in the case of the thermal defect. This leads to 

xl^/3 = x?- + xr = ^. (7.3) 

The one-string defect on the triangular lattice is shown in Fig. ITTT c). It is seen to 
be equivalent to a one-string defect on one hexagonal lattice and no other defects on the 
other hexagonal lattices. The critical exponent is therefore 

xr = 4- = l. (7.4) 
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Figure 17: Defect vortices on the triangular lattice (left) and on the hexagonal lattice 
(right), with the corresponding magnetic charges, (a) Thermal defect, (6) the ±7r/3 bend 
defect, (c) one-string defect, {d)-{f) various examples of two-string defects (see text). On 
the left, the RGB loops on the triangular lattice are shown in thick grey linestyle, and 
the empty bonds by dotted lines. On the right, the string defects on the hexagonal 
lattices (with labels as shown) are given in thin grey linestyle; non-defect directed loops 
are black. The magnetic charges on the two sublattices are denoted by mi and m2. 
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As discussed earlier, the fractal dimension D{ of a loop is such that 2 — Df (the 
codimension) is the critical exponent describing the decay of the probability that two 
distant points belong to the same loop [321 tHH] • In a model of oriented loops, this critical 
exponent can be computed by reversing the direction along one half of the loop, whence 
it becomes identified with X2 as in Eq. ()5.15|) . In the present case, the RGB loops do not 
possess an intrinsic orientation, as discussed above, and the two points must be marked 
in another way. A convenient way is by inserting the bend defect of Fig. m\ b). We 
therefore arrive at 

A = 2 - 4^;/3 = I (7.5) 

This agrees with earlier numerical work jl2]. Note that this is different from Df of the 
n = 1 FPL models on the hexagonal and the square lattices (which both have Df = 7/4). 

There are several different ways of making two strings emanate from nearby, but 
distinct points. The corresponding critical exponents are thus not linked to D^, and 
moreover depend on the exact geometry of the string insertions, even though macroscop- 
ically the two strings appear to have common end-points. We illustrate this comment by 
a few examples. 

Consider first the geometry of Fig. irTT rf) in which the two nearby defects reside 
on the same sublattice and have the same alignment. The defect charge is equivalent 
to that of two nearby one-string defects on one hexagonal lattice, and leads to x*^) = 
x^"^(eo, 2(2B + C)) = 5/4, by Eq. 

A different alignment of the strings is shown in Fig. lT7r e)-( f). If the two defects reside 
on different sublattices, as in Fig. ITTTe). one obtains x^") = x^^^^Bq, (B — C)) = 1/4. If 
they are on the same sublattice, as in Fig. \T7\ f). one has x*"-, = 1/2. 

7.2 The triple RGB model 

A more exotic model can be obtained if we set n = 2, so as to weigh equally the three- 
colouring configurations (rather than the dimer coverings) on each of the three hexagonal 
lattices. Each vertex on the triangular lattice is now visited by three distinct RGB-type 
loops, having respective colours A, B and C . 

A sample configuration is shown in Fig. IT^ a). where the A- loops appear as thick grey 
lines, the 5-loops as thick black lines, and the C-loops as thin black lines. The central 
charge of this model is now c = 3 x 2 = 6. 

Again, we emphasise that the three copies of the n = 2 hexagonal-lattice FPL model 
are completely decoupled, and critical exponents can be obtained along the lines laid out 
above. 

7.3 Staggered models 

Yet other models can be obtained by taking the w ^ oo limit in one or several of 
the superposed hexagonal-lattice FPL models. As we have seen in section 13.51 this will 
eliminate one unit of central charge from each of the concerned models. 
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We recall that infinite w means that the i?o-type edges cannot have the colour A. 
Thus, within the RGB model one obtains critical theories with c = 2, 1 or 0, depending 
on whether the loops are excluded from the Eq edges on one, two or three hexagonal 
lattices. Similarly, within the triple RGB model one gets c = 5, 4 or 3, depending on 
whether the A-coloured RGB loops are excluded from the Eq edges on one, two or three 
hexagonal lattices. In all cases, critical exponents can be worked out as usual. 

8 Conclusion 

In this paper, we have studied the loop model (jl.5|) on the hexagonal lattice subjected to 
a staggered field w. The parameter w breaks down the three-phase coexistence present 
at w = 1 and has an interesting competition with the temperature T. 

We have found that the model has four new critical branches for w 7^ 1, as seen in 
Fig. ^1 Three of these exist for < n < 2 and are denoted by P (Potts model), D' 
(tricritical dense loops) and d' (tricritical dilute loops). Furthermore, there is a fourth 
critical point for n > 2 belonging to the Ising universality class. These are in addition to 
the four old branches known to exist at w = 1: C (compact loops), D (dense loops) and 
d (dilute loops) for \n\ < 2; and the critical three-state Potts model for n > 2. 

It seems that the multicritical branches D' and d' are as generic as the D and d 
branches, since all of these have previously been found on the square and triangular 
lattice 0{n) models. An explanation of this universality is given within the Coulomb 
gas formalism in section 16.31 On the other hand, the point C (compact loops, T = 
and W7 = 1) is specific to the hexagonal lattice. Due to the enhanced symmetry, its 
Coulomb gas is based on a two-component height field, and this is unstable towards both 
the perturbations in T and in w. 

A part of our results can be explained via the embedding of the model ()1.5p into 
the more general triangular-lattice 0(n) model. However, the classification of possible 
types of critical behaviour in this latter model |23j does not appear to be complete. 
In particular, compact loops and the RGB model are interesting special cases of the 
triangular-lattice 0(n) model which are not accounted for in Ref. [22] ■ 

It should be emphasised that even though the 0(n)-type loop models on various 
lattices (hexagonal, square and triangular) share many universal features, they also fail 
to account for some aspects of the low-temperature physics of the original 0(n) spin 
model. For example, the generic symmetry-broken low-temperature Goldstone phase is 
absent [3 CD], and for n = 2 there is no line of critical points below the Kosterlitz-Thouless 
temperature 

Finally, we have also studied constrained compact RGB-type loops on the triangular 
lattice. In particular, we have proved that RGB loops have the fractal dimension = 
3/2. This finding would seem to indicate that these loops are level lines of a Gaussian 
surface. More precisely, we conjecture that the RGB height (defined in section 17. ip 
defines a Gaussian surface. 

Interestingly, our derivation of the value of Df is based on level lines of three different 
four-dimensional surfaces which are manifestly non-Gaussian. It would be interesting if 
one could find a more direct derivation, in terms of the RGB height. 
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